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ABSTRACT 

We use published iJOS'^T observations of the X-ray Nova V1974 Cygni 1992 to test 
a model for interstellar dust, consisting of a mixture of carbonaceous grains and silicate 
grains. The time-dependent X-ray emission from the nova is modelled as the sum of 
emission from a 0-Ne white dwarf plus a thermal plasma, and X-ray scattering is calcu- 
lated for a dust mixture with a realistic size distribution. Model results are compared 
with the scattered X-ray halos measured by ROSAT at 9 different epochs, including 
the early period of rising X-ray emission, the "plateau" phase of steady emission, and 
the X-ray decline at late times. We find that the observed X-ray halos appear to be 
consistent with the halos calculated for the size distribution of Weingartner &: Draine 
which reproduces the Milky Way extinction with Ry = 3.1, provided that the reddening 
to the nova is E{B — V)'^ 0.20, consistent with E{B — V) ^ 0.19 inferred from the 
late-time Balmer decrement. The time delay of the scattered halo relative to the direct 
flux from the nova is clearly detected. 

Models with smoothly-distributed dust give good overall agreement with the ob- 
served scattering halo, but tend to produce somewhat more scattering than observed 
at 200-300", and insufficient scattering at 50-100". While an additional population of 
large grains can increase the scattered intensity at 50-100", this could also be achieved 
by having '^30% of the dust in a cloud at a distance from us equal to '^95% of the 
distance to the nova. Such a model also improves agreement with the data at larger an- 
gles, and illustrates the sensitivity of X-ray scattering halos to the location of the dust. 
The observations therefore do not require a population of micron-sized dust grains. 

Future observations by Chandra and XMM-Newton of X-ray scattering halos around 
extragalactic point sources can provide more stringent tests of interstellar dust models. 



Subject headings: radiative transfer — novae, cataclysmic variables — scattering — 
stars, individual (V1974 Cygni) — dust, extinction — x-rays: ISM 
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1. Introduction 

Interstellar grains scatter X-rays through small scattering angles, and as a result distant X-ray 
point sources appear to be surrounded by a diffuse "halo" of scattered X-rays (Overbeck 1965; 
Martin 1970; Hayakawa 1973). The angular structure and absolute intensity of these scattered 
halos can be measured using imaging X-ray telescopes, thus providing a test for interstellar grain 
models (Catura 1983; Mauche & Gorenstein 1986; Mitsuda et al. 1990; Mathis & Lee 1991; Clark 
et al. 1994; Woo ct al. 1994; Mathis et al. 1995; Predehl & Klose 1996; Smith & Dwek 1998; Witt 
et al. 2001; Smith et al. 2002). Given an accurate model for the dust grain size distribution and 
its scattering properties, observations of scattering halos can also constrain the spatial distribution 
of dust towards a source and the distance to a source, particularly if the emission is time variable 
(Triimper & Schonfelder 1973; Predehl et al. 2000). 

Nova V1974 Cygni 1992, a bright X-ray nova, was observed extensively by the imaging X- 
ray telescope on ROSAT, resulting in the best extant data set for studies of the X-ray scattering 
properties of dust (Krautter et al. 1996). Mathis et al. (1995) compared model calculations to the 
X-ray halo observed 291 days after optical maximum, and argued that the angular structure of the 
observed X-ray halo favored a grain model based on highly porous grains. Smith & Dwek (1998) 
disagreed with this conclusion, arguing that the halo around Nova Cygni 1992 did not require porous 
grains, but was in fact consistent with the scattering expected from a mixture of nonporous silicate 
and carbon grains. More recently, Witt, Smith, & Dwek (2001) reached a different conclusion, 
arguing that the observed X-ray halo around Nova Cygni 1992 requires that the size distribution 
of interstellar dust grains extend to radii a > 2.0//m, with >40% of the dust mass in grains with 
radii a > 0.5/im. 

Weingartner & Draine (2001) and Li & Draine (2001) have recently put forward a physical dust 
model which is in quantitative agreement with the wavelength-dependent extinction of starlight as 
well as the observed spectrum of infrared emission from interstellar dust. The model consists of 
a mixture of carbonaceous grains (including ultrasmall grains with the properties of polycyclic 
aromatic hydrocarbon molecules) and amorphous silicate grains. By appropriate adjustment of the 
size distribution, the model can reproduce the extinction in different regions of the Milky Way, 
and in the Large and Small Magellanic Clouds. Li & Draine (2002) show that the model is also 
consistent with the observed infrared emission from the Small Magellanic Cloud. Here we use the 
observed X-ray halo around Nova Cygni 1992 to test this dust grain model. Since the sightline 
toward Nova Cygni 1992 is presumably typical diffuse interstellar medium, we use the Weingartner 
& Draine (2001, hereafter WDOl) size distribution for Milky Way dust with Ry = 3.1. 

In §3 we review estimates of the distance to and the gas and dust toward Nova Cygni 1992. 
An empirical model for the X-ray emission from the nova is described in §4, with the emission 
modelled as the sum of emission from a hot thermal plasma plus a white dwarf photosphere with 
varying temperature and radius. The methodology for calculation of X-ray scattering by dust is 
presented in §5, including multiple scattering, the effects of time delay, and the calculation of dust 
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scattering cross sections. 

Our results are presented in §6. We test our model using observations from 9 different epochs, at 
radii out to 2000". We show that the WDOl model is in quite good agreement with the observations 
if the dust is assumed to follow an exponential density law and the nova is at a distance of ~2.1 kpc. 
The time delay of the halo relative to the nova is clearly visible at late times when the nova is in 
decline. We discuss the uncertainties associated with possible clumping of the dust into clouds 
along the line-of-sight, and show that agreement with the observed halos can be improved if ~ 30% 
of the dust is concentrated in a cloud ~ 100 pc from the nova. We conclude that the WDOl dust 
model is consistent with the observed X-ray halo toward Nova Cygni 1992, and a population of 
large dust grains is not required. 

The distance estimate to the nova depends somewhat on the assumed density distribution of 
dust, and is considered in more detail in a separate paper (Draine & Tan 2003). 

2. X-Ray Observations 

ROSAT PSPC\m&ges of Nova Cygni 1992 for days^ 258, 261, 291, 434, 511, 612, 624, 635, and 
650 after optical maximum were extracted from the ROSAT aichwc and analysed using the ESAS 
software package (Snowden et al. 1994). This corrects for the effects of the nonuniform detector 
response (particularly the shadowing by struts around 1200") and vignetting, and excludes periods 
of high solar activity, allowing fairly accurate X-ray intensity profiles to be determined out to 
^ 2000 — 3000", particularly when the nova was bright. Prom the observation at day 650, which 
has the deepest exposure and the weakest nova halo, we identified the nine brightest background 
point sources. Excluding these sources,^ the azimuthally-averaged intensities and their statistical 
uncertainties were calculated for all annuli around the nova. 

The diffuse background was estimated for each epoch by averaging over an annulus from 
2800" to 3200", where the angular dependence of intensity was seen to be flat. Backgrounds 
ranged from 5.2 x 10"'^ ct s~^ arcmin"^ (day 261) to 9.9 x 10"'^ ct s~^ arcmin"^ (day 624). This 
scatter is most likely due to uncertainties in estimating intrumental backgrounds in the ESAS data 
reduction process. Background subtraction was therefore done separately for each epoch so that 
these systematic uncertainties woTild have a reduced impact on the determination of the intensity 
of the nova's halo. Note, however, this method of evaluating the background forces the derived 
halos to artificially go to zero at 3000". In Table 1 we list the background-subtracted count rates 
in selected annuli; we do not go beyond 2040" because uncertainties in the background correction 



^Day "258" is the sum of observations on days 255 and 259. Day "650" is the sum of observations on days 647, 
648, 652, and 653. 

^Only IRXS J202742. 6+522920 {6 = 1630") with 0.09 ct/s and IRXS J202742.4-I-523621 (6» = 1510") with 
0.035 ct/s make significant contributions to the halos. 
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dominate at larger angles. The error estimates in Table 1 include only statistical errors plus an 
estimated ±10% error in the background estimate for each epoch. 

Table 1: X-Ray Lightcurve 



epoch 


iV(<50") 


iV(50-100") 


iV(100-300") 


iV(300-960") 


iV(960-2040") 


(day) 


ct/s 


ct/s 


ct/s 


ct/s 


ct/s 


258 


11.47 ± .04 


0.522 ± .009 


0.574 ± .011 


0.93 ± .05 


0.38 ±0.17 


261 


15.10 ±.07 


0.628 ±.015 


0.619 ±.016 


1.06 ± .04 


0.60 ±0.15 


291 


29.42 ±0.11 


1.12 ±.02 


1.23 ± .02 


2.14 ± .07 


1.1 ±0.3 


434 


68.09 ±0.15 


2.26 ± .03 


2.68 ± .03 


4.99 ± .07 


3.0 ±0.2 


511 


66.41 ±0.17 


2.45 ± .03 


2.74 ± .04 


5.14 ± .08 


3.3 ±0.2 


612 


3.68 ± .02 


0.200 ± .005 


0.309 ± .009 


0.81 ± .06 


0.8 ±0.2 


624 


1.37 ±.02 


.093 ± .005 


0.168 ±.010 


0.53 ± .07 


0.6 ±0.3 


635 


0.589 ±.011 


.049 ± .003 


.091 ± .008 


0.31 ± .07 


o.2l°:i 


650 


0.225 ± .004 


.0251 ± .0015 


.041 ± .006 


0.18 ± .06 


o.i3l°o;?J 



Suppose that the point spread function (psf) has a fraction ^{9) of the total counts within an 
angle 0. At the median photon energy E = 480 eV of the detected photons, the ROSAT ^si has 
7(50") ^ 0.9800 and -/(lOO") » 0.9869 (Boese 2000). Suppose that g{e) is the fraction of the halo 
photons interior to 9. Then (neglecting the effect of the psf on the scattered photons) the total 
point source count rate is 

. _ iV(0 - 50") - /3iV(50" - 100") 

" 7(50") - [7(100") - 7(50")] P ' ^ ' 

- aim .2) 

^ ff(100")-5(50") ■ 

A uniform surface-brightness halo would have [3 = 1/3, but halos for continuously-distributed dust 
have g{9) oc 6* for 6* — > (Draine 2003b), corresponding to /3 = 1 if this behavior applied out to 
100". We take /? «i 0.9 as providing a good approximation to the actual models. 

Table 2 presents our derived point source count rate iVps, the estimated halo count rate A^haio 
(including only halo angles < 2040"), and the ratio iVhaio/A'ps of halo counts to point source 
counts. Note that Nhalo/^ps becomes quite large at late times - this is because the point source 
count rate is declining rapidly, but the longer light travel time means that the halo photons were 
emitted from the nova at an earlier (more luminous) time than the unscattered photons. 

At times t < 511, the ratio of halo to point source is approximately constant. Because day 

511 appears to have been preceded by ~100 days of nearly constant X-ray luminosity, the observed 
halo at day 511 can be interpreted as due to steady illumination, and we can estimate the dust 
scattering optical depth to be 

Tsca(^ < 2040") « In [1 + (iVhaio/iVps) J « 0.193 ± 0.004 (3) 
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Table 2: Point Source and Halo Components 



epoch 




Nhalo 




(days) 


ct/s 


ct/s 





258 


11.30 ± .05 


2.6 ± 0.2 


0.23 ± .02 


261 


14.93 ± .08 


3.1 ±0.2 


0.206 ± .013 


291 


29.18 ± 0.12 


5.8 ±0.3 


0.200 ± .011 


434 


67.84 ± 0.16 


13.2 ±0.3 


0.195 ±.005 


511 


65.93 ±0.18 


14.1 ±0.3 


0.213 ± .005 


612 


3.59 ± .02 


2.2 ±0.3 


0.62 ± .08 


624 


1.32 ± .02 


1.4 ±0.4 


1.1 ±0.3 


635 


0.560 ± .012 


0.7 ±0.3 


1.2 ±0.6 


650 


0.208 ± .004 


0.40 ± 0.28 


1.9 ± 1.3 



at a characteristic energy E m 480 eV. Based on the modelling discussed below, we estimate that 
'~91.6=b2% of the halo is at < 2040"; the total X-ray scattering optical depth toward nova Cygni 
1992 is therefore Tsca ~ 0.211 ± .006 



3. Distance and Gas Distribution 

Nova Cygni 1992 (02000 = 20''30'"31.*76, ^2000 = +52<'37'52."9; Austin et al. 1996) has Galactic 

coordinates I = 89.1°, b = 7.82° and is located at a height = Dsinb = 272(1^/2 kpc) pc above 
the Galactic plane. The distance D has been controversial, with recent estimates 3.2 ± 0.5 kpc 
(Paresce et al. 1994), D = 2.1 ± 0.7kpc (Austin et al. 1996), 1.8 ± 0.1 kpc (Chochol et al. 1997), 
and D = 2.6± 0.25 kpc (Balman et al. 1998). 

To calculate X-ray scattering by dust, we require a model for the spatial distribution of dust 
between us and the nova. We take the Sun to be located at the Galactic midplane, z = 0. We 
model the distribution of interstellar gas by an exponential distribution: 

^H"(^^)-"'^ exp(-./>..), = M . (4) 

N^^{D) = Ni^{oo) [1 - exp{-Dsmb/he)] , (5) 

where z is the height above the Galactic plane, and Zi/2 is the height where the density is 50% of 
the maximum value. 

H I 21 cm emission maps of this region (Hartmann &: Burton 1997) indicate a total atomic H 
column density A^(H I) = 2.2 x 10^^ cm^^ within --^5 kpc (integration performed from —50 km s~^ < 
"flsr < +50 km s-i: a weak component at —150km s ""^ < fisr < —50km s ^ is excluded), if the 
emission is optically thin. From the COBE and DIRBE far- infrared maps, Schlegel et al. (1998) 
estimate a total dust column with E{B — V) = 0.412 mag; with the local ratio N}i/E{B — V) = 
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5.8 X 10^^ cm-2 (Bohlin, Savage, & Drake 1978) this corresponds to N^^{oo) = 2.39 x 10^^ cm-^, 
which we adopt as our best estimate in the direction of the nova. 

Based on studies of the vertical distribution of the gas (see Binney Sz Merrifield 1998, Fig. 9.25) 
we adopt a half-density height Zi/2 = 300 pc for the gas at Galactocentric radius Rq < R < 1.03i?o- 
This gives a midplanc density nH(0) = N^^{oo) sinb/he = 0.24cm~^. We will also consider 
Zi/2 = 250 pc and 350 pc for comparison. 

Given an adopted density law, we treat the distance D to the nova as an adjustable parameter 
which determines the column of dust and gas between us and the nova, which in turn determines 
the strength of the scattering halo relative to the X-ray point source. By comparing models with 
different D, we will determine what column density N^^{D) best reproduces the observed strength 
of the scattering halo. 



3.1. Attenuation by Gas and Dust 

Radiation from the nova is attenuated by gas and dust along the line of sight. Prom the ROSAT 
data, Balman et al. (1998) find that the observed flux from the nova is consistent with emission 
from hot plasma plus a white dwarf photosphere, attenuated by absorption by interstellar gas with 
column density A^h ~ 2.1 x 10^^ cm~^ at late times {t > 255 d), but with additional absorption at 
earlier times. 

Our objective is simply to find an empirical description which reproduces the observed count 
rate and energy spectrum of unscattered photons arriving at the Earth. To this end we adopt the 
parameters estimated by Balman et al. The results of Balman et al. appear to be consistent with 
time-dependent absorption by 

(255 d \ ^ 
— pj cm-2 fori< 255d (6) 

= 2.1 X 10^^ cm-2 for t > 255d . (7) 

The decline of with time is presumed to be due to the combined effects of expansion and 
ionization of gas associated with the nova. We will assume that the absorption is due to the 
interstellar contribution [eq. (5)] plus a time-variable 

ATVh (t) = (t) - (8) 

contributed by gas associated with the nova.^ 

We will see below that the observed X-ray halo appears to be consistent with models with 
Njf^'^ « 1.0 X 10^1 cm between us and the nova. The absorption optical depth of the interstellar 



^Note that for all the nova distances and density distributions we consider, < 2.1 x 10^^ cm ^, so that 

ANh > 0. 
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gas is shown in Fig. 1 as a function of energy, for = 1.0 x 10^^ cm~^. Photoelectric absorption 

due to H, He, C, N, O, Ne, Mg, Si, S, and Fe is included, with cross sections calculated following 
Verner and Yakovlev (1995) and Verner et al. (1996), using subroutine phfit2.f written by D.A. 
Verner (1996). 

In the interstellar gas, abundances relative to H arc taken to be 100% of solar for He I, N I, 
Ne I, and S II, 30% of solar for C II, 80% of solar for O I, and 10% of solar for Mg II, Si II, and Fe II 
due to depletion into dust grains. Nova Cygni 1992 did not form dust in the ejecta (Woodward et 
al. 1997). In the circumstellar material we assume no dust grains, and solar abundances for He I, 
C II, N I, O I, Ne I, Mg II, Si II, S II, and Fe II. Solar abundances for He, N, Ne, S are from 
Grcvcssc &; Sauval (1998); for Si and Fe from Asplund (2000); for C from AUende Prieto et al. 
(2002a); for O from Allende Prieto et al. (2002b). 

In addition to gas phase absorption, there is absorption and scattering by interstellar dust 
grains. This is calculated assuming a mixture of carbonaceous grains and silicate grains, with the 
"case A" size distribution found by WDOl for Milky Way dust with Ry = Av/E{B-V) = 3.1, but 
with grain abundance per H nucleon reduced by a factor 0.93, as recommended by Draine (2003a). 
Absorption and scattering cross sections were calculated using dielectric functions which include 
structure near the X-ray absorption edges (Draine 2003b), as described in §5. 

Extinction and scattering optical depths for this dust are shown in Fig. 1 as functions of 
photon energy, for a sig htline with A^^^M ^ i q ^ ;lo21 cm . Below ~ 0.5 keV absorption by the gas 
dominates, since H and He are both assumed to be neutral. At ~ 0.8 keV the dust grains provide 
-^50% of the extinction, and at energies above 1 keV the dust grains dominate the extinction. 

The interstellar matter {N^^) and the circumstellar gas (AA^h) arc taken to have differing 
attenuation properties, since some of the interstellar matter (but none of the circumstellar material) 
is in the form of dust grains. Our estimate for the total attenuation therefore depends on the fraction 
of the total column contributed by the interstellar medium, and this in turn depends on the assumed 
distance D. 

4. X-Ray Spectrum of Nova V1974 Cygni 1992 

Scattering of X-rays by dust grains is a function of the X-ray energy, so it is important to use a 
realistic source spectrum when modelling this phenomenon. Mathis et al. (1995) approximated the 
emission as a blackbody with kT = 22 eV, attenuated by foreground gas with A^h = 4.25 x 10^^ cm~^ 
and solar abundances, while Witt et al. (2001) approximated the spectrum by a delta function with 
E = 400 eV. 

Balman, Krautter &: Ogclman (1998) have recently examined the lightcurve and spectrum of 
Nova Cygni 1992, and find that it appears to be the sum of two separate components: thermal 
emission from the photosphere of the white dwarf, with effective temperature r„d and radius i?wd 



8 




E (keV) 



Fig. 1. — Optical depths due to photoelectric absorption by interstellar gas and extinction by dust, for a sightline 
with iVf ^ = 1.0 X 10^^ cm . Also shown is the contribution to dust extinction due to scattering. Photoelectric 
absorption edges are indicated for C II, N I, O I, and Ne I. 



both varying in time, plus emission from a cooling and expanding thermal plasma. We adopt this 
two-component model for the point-source flux Fy-. 



Fu{t) = te'(t) + F:^'^{t)\ exp [-rabs,gas(i^) - rabs,dust(z^) - r,ca(l^)] 



(9) 



4.1. Thermal Plasma 

The flux Fy^ from the emitting plasma is 

F^' = / n^nndV^^^ (10) 

where Ajy(Tpi), the power radiated per unit volume and unit frequency by an optically-thin thermal 
plasma with kinetic temperature Tp\, is calculated using the Raymond & Smith (1977) code for a 
solar abundance plasma in collisional ionization equilibrium. 
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We take the plasma properties to be given by 

f K 



t < 63d 

10^ K [1 - 0.9(t - 63 d)/42 d] 63 d < t < 105 d 
lO^K t>105d 



■^-jy2 J neTindV = 1.44 x lO^^^picm 
with a correction factor 

A^i = exp 



-5 



0.0925 



1021 cm-2 



(t/94d)2 t<94d 
(i/94d)-2 t>94d 



- 1.0 



(11) 



(12) 



(13) 



depending on the assumed value of N^^. The adopted plasma temperature (eq. 11) is consistent 
with the temperature of this component inferred by Balman et al., and eq. (12) reproduces the 
PSPC count rates at times t < 147 d 



4.2. Nova Photosphere 

Balman et al. found that the emission from the nova photosphere was consistent with the 
spectrum calculated by MacDonald & Vennes (1991) for 0-Ne enhanced white dwarf model atmo- 
spheres. Following Balman et al., we approximate the emission from the nova photosphere using 
the 0-Ne enhanced white dwarf model atmosphere spectra of MacDonald &; Vennes (1991). We 
will assume that the light curve consists of three phases: (1) at times t < tpi, the atmosphere is 
contracting at constant luminosity, with the temperature rising; (2) at times tpi < t < tp2, the 
atmospheric radius remains constant, with the temperature nearly constant (the "plateau" phase); 
(3) at t > tp2 the nova has exhausted its fuel, and the photosphere begins to cool rapidly at constant 
radius. 

The beginning and end of the plateau phase arc not well determined by the ROSAT observa- 
tions. We will assume that the plateau phase began at tpi = 335 d, and ended at tp2 = 565 d. The 
highest observed count rate was aX t = 434 d, at which time we assume an effective temperature 
^^wd = 50 eV; the i^OiS^T point source count rate (see §4.3) on day 434 can be reproduced if 

TT (^^^ = 2.61 X 10-25 sr for t > tpi (14) 



with a correction factor 

A^d = exp 



(15) 



which depends weakly on the adopted value of N^^. With the assumption of i^^d = constant for 
t > tpi, we determine the effective temperature kT^d = 49.84eV on day 511. We assume dT^d/dt 
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to be constant throughout the plateau phase, giving kTpi = 50.16 eV, and kTp2 = 49.72 eV. The 
peak luminosity is reached at tpi: 

L^d = 2.57 X 10^8 (^-^^ ergss-^ . (16) 
At times t < tpi the nova radius is obtained by assuming constant luminosity: 

.(i|l)^2.61xlO-(^)' t„r»<V . (17, 

For t = 258, 261, 612, 624, 635, and 650 days we determine T^d by requiring that we reproduce the 
observed Nps, for constant luminosity at t < tpi, and constant radius for t > tpi. For t < 258 d, we 
assume T (x t. 

Tabulated model atmosphere spectra have been obtained from MacDonald (2002) for temper- 
atures T^d = 1) 2, 3, 3.5, 4, 5, 6 x 10^ K; we estimate the spectrum at intermediate temperatures by 
interpolation. 



4.3. Model Count Rate 

Our aim is to reproduce the ROSAT counts as a function of halo angle. We must first ensure 
that our nova model reproduces the observed count rates for the point-source component. For the 
nova model and interstellar medium parameters described above, we calculate the rate of 0.1 - 2.4 
keV photon detections by the ROSAT PSPC. We use the ROSAT effective area vs. energy from 
Snowden et al. (1994). 

In Figure 2 we show the iJOS'^T point-source count rate calculated for our model for the nova 
and intervening extinction, together with measured count rates. 

In Figure 3 we show the calculated energy spectrum of detected photons at 4 different times. 
We note that the energy spectrum varies considerably over the evolution of the nova. At t = 91 d 
the spectrum is quite hard, being dominated by the thermal emission from hot plasma with kTpi ~ 
3.4 keV. At this time the white dwarf photosphere is relatively cool, and the radiation from it is 
absorbed by intervening H and He. As the white dwarf photosphere contracts and becomes hotter, 
its 0.2-0.7 keV emission comes to dominate the count rate (see, e.g., the spectra for t = 258, 434, 
and 635 d in Figure 3) . 

In Fig. 4 we display the energy spectrum of detected photons at t = 291 d after optical maxi- 
mum, together with the spectra used by Mathis et al. (1995) and Witt et al. (2001). 
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t (days) 

Fig. 2. — Light curve for point source component of Nova Cygni 1992. Solid line: count rate for two component 
model, including attenuation by gas and dust. Broken curves: contributions from thermal plasma and white dwarf 
photosphere. Squares: measured count rates from Krautter et al. (1996) for t < 147 d, and from Table 1 for t > 258 d. 

5. X-Ray Scattering 

The theoretical framework for calculation of scattered halos around X-ray sources has been 
discussed by Mauche & Gorenstein (1986), Mitsuda et al. (1990), Mathis & Lee (1991), and Predehl 
& Klose (1996). The treatment which we present here is fully general, subject to only the following 
approximations: 

1. Polarization is neglected. This is an excellent approximation for X-rays since scattering is 
only appreciable for small scattering angles (< 1°) for which polarization effects are negligible. 

2. Dust grains are assumed to be randomly-oriented. If dust grains are nonspherical and prefer- 
entially aligned, the dust scattering halo would not be azimuthally symmetric. While it may 
be possible to detect this effect in future observational studies, we neglect it in the present 
treatment. Generalization to include this effect would be straightforward but cumbersome, 
and not merited at this time. 

3. We assume that the extinction (scattering by dust plus absorption by gas and dust) along each 
path contributing to the scattered intensity is the same as the extinction along the direct line 
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E (keV) E (keV) 



Fig. 3. — ROSAT PSPC count rates as function of photon energy for two-component model for Nova Cygni 1992, 
including extinction by gas and dust, for = 1.15 x 10^^ cm^'^. 

from observer to source. Since the scattering halos are small (~90% of the scattered photons 
within ~ 0.5°) this is quite a good approximation. Towards Nova Cygni 1992 there are 
variations in A^^^'^(oo) of ~ 20% on 0.5° scales, particularly in a direction perpendicular to 
the Galactic plane (Hartmann & Burton 1997). In the optically-thin limit (i.e. for hiy > 1 keV, 
see Fig. 1), such a density gradient has no effect on the azimuthally-averaged halo intensity, 
relative to uniform extinction with a column equal to the average over the halo scale in the 
inhomogeneous case. When the optical depth is ~ 1 (i.e. for hi^ ~ 500 eV), the differential 
absorption makes a small (~ 1%) effect on the azimuthally-averaged profiles compared to the 
uniform case. 

4. It is assumed that photons scattered by more than 90° may be neglected: we only consider 
scattered photons travelling away from the source (the "outward-only approximation"). 

5. The dust is assumed to be distributed spherically-symmetrically around the source. For 
example, there could be a spherical cavity surrounding the source.^ Outside of this cavity, 
when we use a prescription (e.g., exponential density law) for the dust density, we take this 



^In our calculations we assume a small cavity around the source to avoid numerical difflculties with divergence in 
the integrand of eq. (19) for ^ and y ^ 1. 



13 



102 



m 




m 
-i-j 


IQi 


a 




ID 




O 




o 








H 


1 


T3 


\ 








W 


0.1 



10- 



0.1 



Nova Cygni 1992 
day 291 




WSDOl 
MCFK95' 



0.2 0.3 0.5 0.7 
E (keV) 



Fig. 4. — ROSAT PSPC count rate spectrum as function of photon energy for day 291. Present model (solid 
histogram) is compared with count rate spectrum adopted by Mathis et al. (1995) (MCFK95) and the representative 
photon energy used by Witt et al. (2001) (WSDOl). The drop in count rate at 284 eV is due to absorption by C in 
the 7? 05^ T window. 



to be the density along the sighthne to the nova; the assumed spherical symmetry provides 
densities away from the sightline. Within the small halo angles < 1° of interest, the resulting 
densities are close to what would have been obtained for the plane-parallel density structure 
which is the basis for our original density estimate. The assumed spherical symmetry implies 
azimuthal symmetry around the sightline to the source. Predehl & Schmitt (1995) found that 
X-ray halos in i? 05"^ T observations did not show any detectable azimuthal asymmetries. 



Our formulae are presented without the small-angle approximations which are frequently found in 
discussions of X-ray scattering by dust. 

Consider a point source of specific luminosity Ly at a distance D from the observer, and let 
r be the distance along the line from the observer to the source. We assume the the dust size 
distribution and composition are the same everywhere, but the dust spatial density p may vary 
along the line of sight. 

Let Tabs and Tgca be the total optical depths for scattering and extinction from source to 
observer. At a distance r = xD along the line from observer to the point source, and at "retarded" 
time t (time measured relative to the arrival of a fiducial pulse emitted by the source) , the intensity 
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Fig. 5. — Geometry [upper panel] for calculation of normalized single scattering intensity h{x, 6) and [lower panel] 
for calculation of multiple scattering intensity In{x,6) given 7„_i(j/,a) for y > x. The observer is at a; = 0. 

of n— times scattered photons arriving from an angle 6 is given by 

InAr,t,e) = ^e--'^^^ In{x = r/D,t,e) . (18) 

Eq. (18) serves to define In{x,t,9), which depends on the dust distribution along the line-of-sight, 
but not on the quantity of dust present. The normalization is such that in the approximation of 
small scattering angles, for a steady source one has J ln{0,t, 6)2Tr sinOdO = 1. 



5.1. Single Scattering 

For single-scattering (see Fig. 5) the dimensionless intensity function Ii is given by 

' ' ^ cos^y^ ^ (1 - y)2 + (y - a;)2 tan2 e L^{t) ^ ' 

where 

y' = l- Y^(l-y)2 + y2tan2^ (20) 

and the dimensionless density is 

~p{x)^ "^'Z"''^ . (21) 
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The dimensionless differential scattering cross section is 

a{es) = — , 

^sca ' 

where da/dO, is the differential scattering cross section for scattering angle 

Qs{x,y,9) = i;o{x,y,9) 

where 



il^o{x, y,6) = + arctan 



{y — x) tan 9 



1 - 



and the time delay 



Sh{x,y,e) = ^l^ + {l-y) 



C COS ( 



tan^e 



1/2 



-(1-x) 



D {1 — x){y — x) 



c (l-y) 2 



28.0 d 



D \ (1 - .x)(y-.x) / 6 



For a steady source, 



dy 



2.0kpcy (l-y) VlOOO" 



p(2/')^(0. 



cos 9 '"^ (1 — y)2 + (y — a;)^ tan^ ' 



(22) 

(23) 
(24) 

(25) 
(26) 



(27) 



5.2. Multiple Scattering 

For a steady source, the intensity of multiply-scattered photons is obtained from the recursion 
formula (Mathis & Lee 1991; Predehl & Klose 1996) 

Il{x.O) = 2 dyp{y') da sm a r^_,{y, a) d^ a{Qs) (« > 2), (28) 

cos t) Jo Jq 

where the scattering angle 

05(x, y, 0, a, /3) = arccos (cos a cos 'i/'o + sina cos/3sin^o) • (29) 

The upper limit 7r/2 on the integral over a corresponds to the "outward-only" approximation 
mentioned above. For a nonsteady source, it is possible to allow exactly for the time delays on 
different light paths, but the formalism becomes cumbersome and the computations burdensome. In 
the present application multiple scattering is a relatively minor effect, and time delay is a relatively 
minor correction, so for n > 2 we make only an approximate correction for time delay by assuming 
that the first n — 1 scatterings occur exactly midway between the source and the last scattering: 

rd(3a{Qs) , (30) 
Jo 



in{x,t,9) 



n 



cos 9 



dy p{y') da sin a i^_i{y, a) 
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6t2{x, 9, a) 





X 




cos 6 


D 


X 



+ 



[(1 - x)"^ + x2 tan2 9] 



1/2 



2c 



(1 



-9^ + (1 



cos a 
x)a'^ 



(31) 
(32) 



To calculate the intensity l„(0,t, 0) of n-times scattered photons at a; = 0, we precalculate 
and tabulate I^_i{y,9) at selected values of y and 9, and then obtain by interpolation when 
evaluating eq. (30) for /„(0, x,t). 
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Fig. 6. — WDOl size distributions for carbonaceous grains (including PAHs) and amorphous silicate grains for 
Milky Way dust with Ry = 3.1 and C/H=60ppm in PAHs. Total volume of carbonaceous grains and silicate grains 
per H is 2.26 x 10"^^ and 3.96 x 10"^'^ cm^, respectively. 



5.3. Dust Model and Scattering Cross Sections 

For the line of sight to Nova Cygni 1992, we assume that the dust is "average" diffuse cloud dust, 
and we adopt the dust model developed by WDOl for Ry = Ay / E{B — V) « 3.1 and C/H = 60 ppm 
in polycyclic aromatic hydrocarbons (PAHs), except that, following Draine (2003a), the abundances 
of all grain components (relative to H) are reduced, by a factor 0.93; with this adjustment, the 
grain model reproduces the estimated extinction per H nucleon. The model consists of a mixture 
of carbonaceous grains and amorphous silicate grains. As discussed by Li &: Draine (2001), the 
carbonaceous grains have the properties of PAH molecules when they contain < 10^ C atoms, 
and the properties of graphite particles when they contain > 10^ C atoms. By altering the size 
distributions of the carbonaceous and silicate particles, the dust model appears to be able to 
reproduce observed extinction curves in various Galactic regions, and in the Large and Small 
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Magellanic Clouds. This dust model, when illuminated by starlight, produces infrared emission 
consistent with the observed emission spectrum of the interstellar medium (Li & Draine 2001, 2002). 
The size distribution for this dust model is shown in Fig. 6. For assumed densities of 2.2gcm~^ 
for the carbonaceous grains and 3.8gcm~^ for the MgFeSi04 silicate grains, this corresponds to a 
dust-to- gas mass ratio of 0.008 

The scattering properties of a spherical target of radius a are determined by the complex 
refractive index m(A) and the dimensionless size parameter x = 27ra/A. Many early papers on 
X-ray scattering by dust employed the "Rayleigh-Gans" approximation to calculate the differential 
scattering cross sections. Smith &: Dwek (1998) showed, however, that the Rayleigh-Gans validity 
criterion \m — l\x <^ 1 is not satisfied for interstellar grains at energies <1 keV. In the present work 
we assume spherical, homogeneous grains and use exact Mie scattering theory (see, e.g. Bohren & 
Huffman 1984) to calculate differential scattering cross sections da/dQ, employing a computer code 
based upon the subroutine MIEVO written by Wiscombe (1980, 1996), modified to use IEEE 64 bit 
arithmetic. Wiscombe's code is accurate for "size parameters" x = 7800(a///m)(keV//ii^) < 2x 10"^. 
For X > 2x 10^ we calculate da/dO. using "anomalous diffraction theory" (van de Hulst 1957), since 
the validity criteria x ^ 1 and |m — 1| ^ 1 are satisfied. The accuracy of the scattering calculations 
is illustrated in Fig. 7, where we show the differential scattering cross section calculated for a case 
with X = 2 X lO'' - the Mie theory calculation and the anomalous diffraction theory calculation 
are indistinguishable. Since the two calculations are entirely different in approach, the agreement 
confirms that both are accurate. 

For the amorphous silicate grains and the carbonaceous grains we employed the dielectric 
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functions estimated by Draine (2003b) for olivine MgFeSi04 and graphite, respectively. Figure 

8 shows the differential scattering cross section per H nucleon for the dust mixture of WDOl for 
Milky Way dust with Ry = 3.1 (with grain abundances reduced by a factor 0.93). For this mixture 
we see that the scattering for scattering angles 0^ < 1000" is dominated by grains with radii in 
the 0.1 — 0.4/im range - larger grains contribute only ~20% of the total scattering cross section 
for @s ^ 200", and < 1% for Qg ^ 500". We also see that silicate grains provide ~60% of the 
scattering for 6^ < 300", increasing to ~80% at 6^ > 1000". 

It is well-known that the extinction curve varies from one region to another, which is presumed 
to be due to changes in the size distribution of the dust. In addition to the size distribution which 
reproduces the standard Milky Way diffuse cloud extinction curve with Ry = 3.1, WDOl have 
constructed a size distribution consistent with an extinction curve with Ry = 5.5. Such flat 
extinction curves are found in dense regions, and require a shift of the dust size distribution toward 
larger sizes. The X-ray scattering properties for dust with Ry = 3.1 and 5.5 are compared in Fig. 

9 - we see that Ry = 5.5 dust is slightly more forward-scattering than Ry = 3.1 dust. The cross 
section for a scattering angle of < 300" is larger by about a factor 1.4. 

The scattering also depends on the X-ray energy. In Fig. 9 we show the differential scattering 
properties of Ry = 3.1 dust at energies ranging from 0.1 to 1 keV. As the energy increases, the 
dust becomes more forward-scattering. 
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Fig. 8. — Heavy solid curve: differential scattering cross section per ¥1 at hv — 450 eV for the WDOl model for 
Milky Way dust with Rv = 3.1 and C/H=60 ppm in PAHs. Also shown are separate contributions of silicate and 
carbonaceous grains, and contributions from different size ranges. 
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Fig. 9. — Upper panel: differential scattering cross section at _E = 450 eV for WDOl model for Milky Way dust 
with Rv = 3.1 (diffuse cloud average) and Rv = 5.5 (dense cloud). Lower panel: differential scattering cross sections 
at 8 energies for WDOl model for Milky Way dust with Rv = 3.1. 



6. Results 
6.1. Calculations 

Using the method described in §5, the dust-scattered halo was calculated at each of 9 epochs, 

and the energy-dependent ROSAT PSPC psf of Boese (2001) was used to model the contribution 
of the unscattered photons to the image. We include the contributions to the scattered halo from 
singly and doubly-scattered photons. For the scattering optical depth Tgca < 0.3 for E > 250 eV 
(see Fig. 1), photons scattered 3 or more times contribute only a fraction 

3! 'sea ^ 4! ^sca ^ - ^ 

gTsca — 1 ^ ' 



of the total halo counts. 
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Fig. 10. — X-ray intensity profile from Nova Cygni 1992, 434 days after optical maximum. Solid histogram: raw 
ROSAT PSPC data, in 10"bins. Triangles: intensities after processing by ESAS software (see text). The background 
is determined by averaging the intensity in the annulus from 2800-3200". Squares: nova intensity after background 
subtraction. The error bar around each point shows the 3 — a statistical error but does not include any uncertainty due 
to background subtraction. Solid curve: model intensity profile [psf (dotted line) plus dust-scattered halo (dashed 
line)] for nova at D = 2.1 kpc. The model is in generally good agreement over the full 40-3000"range where the 
scattered halo can be measured, but overpredicts the scattered intensity by ~40% at 6* ~ 250", and underpredicts 
the observed intensity at 6 < 100" . 

We use 71 energy bins extending from 0.1 to 2.4 keV, chosen so that absorption edge structure is 
weh-defined. The dust size distribution is treated using 50 size bins spanning the range .00035/um < 
a < 1.0/im (except for the "big grain" case discussed below, where we use 55 bins running from 
.00035/im < a < 2.0/im). We calculate the scattering halo at 80 different halo angles < < 
lO'^arcsec. The dust scattering properties i5"(6s) are precalculated at 115 different scattering angles 
Qg, with interpolation used during evaluation of Ii [eq. (19)] and I2 [eq. (30)]. 

For computation of the second-order scattering, Il{x,9) is evaluated at 40 different locations 
X, with interpolation used for evaluation of l2- We have verified the accuracy of the numerical 
quadratures by doubling the number of spatial steps and confirming that the results do not change 
significantly. 

Our preferred exponential model, with an assumed distance D = 2.1 kpc, is shown in Figure 10, 
including the separate contributions of the point source (broadened by the psf) and the scattered 
halo. 
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Fig. 11. — Multi-epoch X-ray intensity profiles from Nova Cygni 1992. Squares: observed intensity profile; the error 
bar shows the 3 — cr statistical uncertainty. Solid lines: models consisting of psf plus dust-scattered halo for exponential 
density law with zi/2 = 300 pc and nova at an assumed distance of D = 2.1 kpc. This distance, corresponding to 
Nli^ = 1.15 X 10^^ cm ^, or E{B — V) — 0.20 mag, gives the best overall match to the data. Dotted lines: the same 
models but for D — 1.0 (lower line) and 3.0 kpc (upper line) (corresponding to = 0.64 and 1.46 x 10^^ cm~^, or 

E{B -V)^ 0.11 and 0.25 mag. 

6.2. Comparison with Observational Data 

The radial distributions of ROSAT PSPC counts are shown in Figure 11 for the nine epochs 
for which high quality i^OS'^T data are available (see §2). 
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As discussed in §2, our method of evaluating the background forces the derived halos to ar- 
tificially go to zero at 3000". Thus the observed halo intensities arc only reliable for intensities 
significantly greater than the true halo intensity at 3000" and for this reason we restrict the useful 
comparisons of models with data to the region inside ~ 2000". An example of the raw, processed, 
and background-subtracted data is shown in Figure 10 for day 434, which had the highest count 
rate (67.8 cts~^) and number of counts (~200,000') of unscattered photons from the nova. 

We also show models calculated for several values of the nova distance D, ranging from D = 1.0 
to 3.0 kpc (holding fixed the assumed exponential gas distribution, with 2:1/2 = 300 pc). As D is 
increased, the dust scattering optical depth Tgca to the nova increases, leading to an increase in the 
intensity of the scattered X-ray halo. Note that as D is varied the observed point source flux is 
held constant by scaling the intrinsic nova luminosity. 

To quantify the goodncss-of-fit, for each of the epochs k = 1-9 we divide the range 50-2040" 
into annuh, j = 1-36, with widths M = 10" for 50-140", 20" for 140-360", 60" for 360-960", and 
180" for 960— 2040" (the nova data are plotted in these intervals in Figures 10 and 11). Let if^ be 
the observed intensity (counts s~^ arcmin"^) including the background for annulus j and epoch 
k, let ff^ be the estimated background intensity in annulus j at epoch k, and let I^S^ be the 
intensity calculated for the model at epoch k. We use a goodness-of-fit metric 

9 36 [robs jbkg rmod 

k=\ j=i ^j'k 

al, ^ lf,V(f^. Ai,) + 0.01 [l^l^) ' + 0.04 - Jbkg) ' + 0.04 [if^^) ' . (35) 

The first term in eq. (35) is due to the statistical uncertainty in the number of photons counted (J7j 
is the solid angle of annulus j, and At^ is the exposure time for epoch k). The second term allows 
for an estimated ±10% uncertainty in the estimated background level. The third and fourth terms 
are somewhat arbitrary, but are introduced to avoid overly weighting the regions of the halo where 
the photon statistics may be good, but where there may still be unknown systematic errors in the 
observations, as well as systematic errors in the models due to inaccuracies in the adopted nova 
light curve, nova spectrum, grain dielectric functions, etc. Note, for example, that even though 
the count rate appears to be systematically rising between days 255 and 292 (Krautter et al. 1996) 
the measured count rate on day 292 is 11% smaller than the count rate measured on day 291, 
indicating either a variation in instrumental response or complexities in the nova emission which 
are not allowed for in our model. With defined by eq. (35) , a single data point cannot contribute 
more than 25 to x^. 

Fig. 12 shows ^ ^ function of for three exponential density laws. For zi/2 = 300 pc, 
is minimized for column density pa 1.15 x 10^^ cm~^, corresponding to a distance D 2100 pc. 

This model is shown in Figs. 10 and 11 by the solid line. At shorter distances is smaller and 

the model halo is systematically too weak, while for larger distances is larger and the model 
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halo is sytematically too strong. The best-fit model still shows systematic discrepancies with the 
observations - it can be seen from Fig. 11 that the model tends to overpredict the halo intensity 
for 200-500" for t < 511 d, and underpredicts the halo intensity for 50-100" at all times - but the 
model halo profiles are in generally good agreement (± < 40%) with the data at all epochs. 

Fig. 12 shows that exponential density laws with different values of 2^/2 nearly degenerate: 
the best-fit distance in each case corresponds to the distance at which the column density N}i{D) ~ 
1.15 X 10^^ cm~^, OT D ^ 2100(2;i/2/300pc) pc. For a steady source these solutions would in fact be 
perfectly degenerate - the difference between models with the same A^^^^ but different Zi/2 is due 
to the variability of the source and the different time delays for models with different D. The 
curves for the 3 different exponential distributions do differ slightly, but evidently the differences in 
time delay have a relatively small effect on the overall x^- This is because the effect of time delays 
on the halos is only prominent at late times when the lightcurve is declining and at large angles 
where the fluxes are close to the level of the background (see Fig. 13). 

Fig. 13 shows X-ray halos at each epoch for our standard model, as well as the results for the 
same model but with time delays set to zero (e.g., infinite speed of light). Time delay effects are 
minimal at days 434 and 511 because the light curve (Fig. 2) is relatively constant for >50 days 
prior to those two epochs; at earlier epochs {t < 291 d) the effect of time delay is to reduce the 
intensity at large angles, and at t > 612 d the effect of time delay is to increase the intensity at 
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Fig. 13. — Effect of time delays and of adding a cloud close to the nova: Same as Fig. 11 but showing D = 2.1 kpc 
model (dotted line) with same model calculated neglecting time delay effects (dashed line). At early times {t < 291 d) 
the effect of time delay is to lower the intensity at large angles; at late times (lower panel) the effect of time delay 
is to increase the halo intensity at large angles. A "cloud + exponential" model, also for D = 2.1 kpc, with 70% of 
dust in the exponential distribution and 30% in a cloud at y = r/D = 0.95 is shown by the solid line (see also Figs. 
14 & 15 and §6.4) 



large angles. 
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6.3. Reddening 

Our modelling of the X-ray lialos favors a column density N^^"^ ~ 1.15 x 10^^ cm~^ in order 
to get the overall halo intensity correct; this column of gas and dust corresponds to a reddening 
E{B — F) PS 0.20 mag. How does this compare to other reddening estimates? 

Bargcr et al. (1993) measured the Hcc/H/? intensity ratio at different times. The ratio declined 
with time, presumably due to the nebula becoming optically thin in the Balmer lines. The last 
measurement reported by Barger et al. {t = 450 d after optical maximum), was /(Hq;)/J(H^) = 
3.47; Mathis et al. (1995) make a reasonable extrapolation to an asymptotic value 3.31, from which 
they estimate E{B — V)= 0.19 mag. If the nebula is optically thick in the Balmer lines, radiative 
transfer effects can lead to an increase in the Ha/H/J ratio relative to the "case B recombination" 
value which is assumed; this would cause E{B — V) to be overestimated. Similarly, at high densities 
there can he n = 2 ^ 3 coUisional excitation, which would again increase Ha emission and lead 
to overestimation of E{B — V). Thus the Hq/H/3 estimate of E{B — V) should be regarded as an 
upper bound. Since the Ha/H/3 estimate of E(B — V) is in agreement with our estimate, it appears 
that radiative trapping effects do not appreciably affect the Ha/H/3 ratio at t > 450 d. 

Austin et al. (1996) estimate E{B -V) = 0.40 ±0.07 using the He H 1640/4686 line ratio, and 
0.38 ± 0.07 using the [Ne IV] 1602/4724 line ratio. In principle these line ratios should be reliable 
reddening indicators.^ However, the He II and [Ne IV] reddening estimates both rely on lUE 
photometry in the 1600-1640A range, and we note that the He II 1640/4686 and [Ne IV] 1602/4724 
ratios measured by Austin et al. at 5 epochs are strongly correlated, which may indicate calibration 
errors in the lUE spectrophotometry. 

If the interstellar gas column on the sightline to the nova substantially exceeds 1.15 x 10^^ cm~^ 
and the true reddening significantly exceeds 0.20 mag, the WDOl dust model will be disfavored: if 
the column density of interstellar dust and gas is appreciably increased, the predicted halo intensities 
become too large at most angles and most epochs, and the goodness-of-fit suffers (sec Fig. 12). We 
conclude that the WDOl dust model is incompatible with E(B — V)> 0.3 in smoothly-distributed 
dust toward the nova. If E{B — V) is indeed as large as 0.36 (the value recommended by Austin 
et al.) then either an appreciable fraction of the reddening must be contributed by dust which is 
located close to the nova (e.g., dust associated with the nova, or interstellar dust within ~50 pc of 
the nova) which would contribute only to very small halo angles, and would be indistinguishable 
from the instrumental point spread function) or the WDOl model must be rejected. However, it 
seems most likely that the reddening is close to E{B — V)^ 0.19 mag, the value obtained from the 
Ha/H/3 ratio.6 



®The [Ne IV] line ratio should be essentially independent of nebular conditions (Draine & Bahcall 1981). If Hq/H/3 
is not appreciably affected by radiative trapping effects, then the He II 1640/4686 ratio would be expected to also be 
consistent with case B recombination theory. 

®From the similar colors of nova V1974 Cygni 1992 and nova V382 Velorum 1999 (for which the H column has 
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6.4. Dust Distribution Along the Line-of-Sight 

The profile of the scattering halo depends on the distribution of dust along the line of sight 
(e.g., Predehl &; Klose 1996). So far we have investigated smooth exponential models. Other 
smooth distributions, such as gaussian or sech^(z) models, do not change the halos significantly, as 
long as the intervening H I column remains close to 1.15 x 10^^ cm~^ (this might require the nova 
to be at a different distance). However, since the total reddening to the nova appears to be just 
E{B — V) ^ 0.19 mag, it would not be implausible for much of this to be contributed mainly by 
one or two diffuse clouds. In Fig. 14 we show models where the same column density of dust is 
located in a single cloud located at y = 0.10, 0.25, 0.5, 0.75, or 0.90 of the distance to the nova. 
We see that the halo profile is quite sensitive to the location of the dust along the line of sight: at 
6 = 100", for example, the halo intensity in the y = 0.10 case is only 6% of the intensity for the 
y = 0.90 case. For a fixed halo angle, there are three separate effects as the dust cloud is moved 
from small y to large y: 

1. The inverse square law causes the intensity of the radiation illuminating the dust grains to 
increase. This acts to increase the surface brightness of the dust cloud. 

2. The required scattering angle 6s increases [see eq. (23) and (24)]. Since da/dO, tends to 
decrease for increasing 0^, this acts to reduce the scattered intensity. 

3. The time delay 5ti increases as y increases. For y = 0.25, 6ti = 9.4d(D/2.0kpc)(^/1000")^, 
while for y = 0.75, 6ti = 84.0d(D/2.0kpc)(^/1000")2. At times when the nova light curve is 
rising (e.g., day 291) this leads to weakening of the halo at large angles ^ as y is increased. 

At small halo angles, the near-constancy of da/dO, at small Qg (see Fig. 8) causes the first effect to 
dominate, so the halo intensity increases as x increases (e.g., at 100", the halo is ~5 times brighter 
for y = 0.75 vs. 0.25). At large halo angles, the rapid decline in da/dQ at large Qg also becomes 
important, as does the time delay if the lightcurve is rising or falling. 

In Fig. 15 we show a model where the total column density to the nova is held at 1.15 x 
10^^ cm~^, but 30% of this is concentrated in a single "cloud" at y = 0.95, with the remaining 
70% of the dust in an exponential distribution with zi/2 = 300 pc. The agreement between model 
and observation is excellent. We see that the dust in the cloud near the nova raises the halo 
intensity at ^ < 150", in accord with observations, and removing 30% of the dust from the general 
exponential distribTition acts to lower the intensity at 200-500", where the model was previously 
somewhat above the observed halo. We show the above "cloud + exponential" model for all 
epochs in Fig. 13. At each time, agreement with the data is improved relative to the best smooth 
model, and is in general excellent. Note that this model is not necessarily the optimum "cloud + 



been measured to be Nh « 1.2 x 10^^ cm ^) Shore et al. (2003) conclude that the reddening to nova V1974 Cygni 
1992 is probably in the range 0.2-0.3. 
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Fig. 14. — Comparison of the X-ray halo produced by smooth exponentiaUy-distributed dust with z^ii ~ 300 pc 
and the nova at D = 2.1 kpc, with models with the same nova distance D and column density A'h but in a single 
cloud a.ty = r/D = 0.1, 0.25, 0.5, 0.75, and 0.90. 

exponential" model. In addition to raising the inner and lowering the outer halo, the dumpiness in 
the distribution imprints structure into the late-time (declining light curve) halos at larger angles 
(see Fig. 13). 

Unfortunately, we have little way of determining the actual distribution of the gas towards the 
nova, since the nova's location at I = 89.1° renders radial velocities unusable for determinations of 
distances within ~ 1 or 2 kpc. In view of these uncertainties, we conclude that the observed X-ray 
scattering halo appears to be consistent, within the uncertainties, with the WDOl dust model. 

6.5. Very Large Grains? 

From Fig. 8 we see that for the WDOl dust model, the X-ray scattering is dominated by the 
grains with radii 0.1/im ^ a < 0.4/im. Grains larger than 0.4/xm contribute less than 20% of 
the scattering even at small scattering angles Qs ^ 100", and less than 1% at scattering angles 

> 1000". 

A recent paper by Witt et al. (2001) has argued that the observed X-ray halo around Nova 
Cygni 1992 requires that the dust grain size distribution include significant mass in large dust 
grains. Witt et al. calculated the dust scattering halo at t = 291 d for a monochromatic /iz^ = 450 eV 
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Fig. 15. — Halo calculated for nova a.t D = 2.1 kpc and (1) standard dust in an exponential distribution (solid line) 
(2) standard dust plus additional equal mass in large grains with a — 2.0/im in an exponential distribution (dashed 
line) (3) standard dust, with 70% of dust in an exponential distribution and 30% in a cloud ai y = r/D — 0.95. 

steady point source. They assumed A^h = 2.1 x 10^^ cm~^, and found that a conventional "MRN" 
{dn/da oc a~^-^) size distribution for a < 0.25/xm (Mathis, Rumpl, Sz Nordsieck 1977), with the dust 
distributed uniformly along the line-of-sight, reproduced the observed halo intensity at 9 = 100", 
but overpredicted the halo intensity by a factor ~2 for 300 - 800"; Witt et al. did not consider halo 
angles larger than 800". In order to suppress the scattering at large angles, while maintaining the 
observed halo intensity at small angles, they favored size distributions with the dust mass shifted 
into larger grains. However, the modifications to the grain size distribution proposed by Witt et 
al. are inconsistent with the average optical and ultraviolet extinction curve for the Milky Way/ 
and there is no reason to assume that this particular sightline has an anomalous extinction law. It 
is our contention that the excess scattering at 300-800" found by Witt et al. is due primarily to 
assuming too much dust on the line-of-sight - our preferred models have Nn ^ 1.15 x 10^^ cm~^, 
only ~ 55% of the 2.1 x 10^^ cm^^ column assumed by Witt et al. 

We have taken our best-fitting smooth exponential model - which uses the standard WDOl 
size distribution which reproduces the standard Ry = 3.1 extinction curve - and modified the grain 
population by adding very large silicate and carbonaceous dust grains. These grains are arbitrarily 



^Witt et al. note that this size distribution has Rv = Ay /{Ab — Ay) ~ 6.2, exceeding the largest values observed 
in diffuse regions. Dust in diffuse regions typically has Rv = 3.1 ± 0.5 . 
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given a radius a = 2/xm, and abundances such that the total sihcate mass and total carbonaceous 
mass is doubled: this trial dust model has 50% of the dust mass in a = 2.0/im grains. Of course 
this dust model now has twice as much dust mass per unit area toward the nova as the original 
model. 

The results of this model are shown in Fig. 15. The additional dust grains do increase the 
strength of the scattered halo, but only slightly. The effect is most noticeable in the range 50- 
100", where the halo intensity is increased by up to 20%. We see, however, that even with this 
unreasonably large mass in large grains, this model is not superior to the model discussed in §6.4 
where 30% of the dust is assumed to be located in a cloud at a distance ~ 100 pc from the nova. 
Adding large grains has negligible effect at large angle halos, whereas placing a cloud close to the 
nova reduces the intensity at large halo angles (since the amount of dust at large distances from the 
nova is reduced); this reduction appears to improve agreement with the observations. We conclude 
that while we cannot rule out a population of large grains, these grains would have only a modest 
effect on the scattered halo, and the observations do not require their presence. A more important 
and plausible effect on the halo comes from the presence of structure in the ISM. 

7. Discussion 

Figure 16 shows the importance of various features of the modelling on the intensity of the 
X-ray halo. The curve labelled "ref" is the single-scattered intensity for a model where the source 
is radiating steadily at a single energy hv = 400 eV, with uniform dust density between observer 
and source. Other curves show the effect of including multiple scattering, time delay, and a realistic 
nova spectrum (same point source count rate) , and replacing the uniform dust distribution with an 
exponential distribution. 

• Doubly-scattered photons add ~10% to the intensity at ^ 600", and ~20% at 6* > 1200". 

• Changing from a uniform to an exponential dust distribution reduces the intensity for 6 < 450" 
(less dust close to the source), and increases the intensity for 9 > 450" (more dust far from 
the source). The increase is ~10% for 6 > 800". 

• Replacing the monochromatic spectrum with a realistic spectrum reduces the flux at large 
angles; the adopted hv = 400 eV is probably slightly low compared to the more realistic 
spectrum (see Fig. 4), resulting in slightly more large angle scattering than for the more 
realistic spectrum. 

• Since the light curve is rising at day 291, including time delay leads to a reduction in the 
intensity, particularly for larger halo angles for which the time delays are larger. At 1000", 
inclusion of time delay reduces the intensity by ~30%. 
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Some of these effects increase, others decrease, the intensity. When all are included together, the 
intensity for day 291 is reduced at angles > 60", with a reduction by ~30% at 1000". 




6 (arcsec) 

Fig. 16. — TvO^I for reference model: hv — 400 eV, uniformly-distributed dust, single scattering only, time delay 
neglected (curve labelled "ref"). Other curves show effect of adding multiple scattering, changing from monochromatic 
to realistic nova spectrum; including effects of time delay for scattered photons, for D = 2.1kpc; and changing from 
a uniform density to exponential density profile with 21/2 = 300pc. Curve labelled "all" is exponential density model 
with full spectrum, multiple scattering, and time delay. 

For smoothly-distributed dust and a standard dust-to-gas ratio, our modelling strongly favors 
a column density ~ 1.15 x 10^^ cm~^ between us and the nova; this is only ~50% of the 

estimated total gas column in this direction. This means that the nova must be close enough to 
have ~50% of the column density beyond it. For the exponential distribution (eq. 4), this places 
it at a distance ~ 2.1(2x/2/300pc) kpc. Since the actual variation of gas density as a function of 
height is not well known, and since much of the gas and dust is likely to be contributed by discrete 
clouds, it is possible for the nova distance to be as large as 2.5 kpc and still have ~50% of the gas 
and dust beyond the nova. 

Most estimates of the reddening to the nova have been larger than the value E{B — V) = 
0.20 mag which we favor. Austin et al. took the weighted mean of 5 different methods, and obtained 
E{B—V) = 0.36ib0.04 mag. Our model with E[B — V) = 0.20 mag has good overall agreement with 
the observed halo intensities, but an increase in E{B — V) to 0.36 would result in halo intensities 
almost twice as strong as observed. If the reddening were shown to be E{B — V^) ^ 0.3 toward Nova 
Cygni 1992, this would rule out the WDOl dust grain model which we have used here i/the dust is 
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smoothly distributed. However, we note that the reddening estimated from the Balmer decrement 
at late times is consistent with our estimate. 

Furthermore, we show that a plausible modification of the spatial distribution of the dust can 
produce good agreement between the observed and calculated halos: if ~ 30% of the dust (i.e., 
E{B — V) = 0.06) is located in a "cloud" within ~ 100 pc of the nova, the calculated scattering 
halo is in good agreement with observation. Since such a spatial distribution is not improbable, 
the observed X-ray halo around Nova Cygni 1992 does not require the existence of a population of 
large grains. Recent observations by Chandra of the scattered halo around GX 13+1 also do not 
support an additional population of large grains (Smith et al. 2002) . 

While we have shown that the WDOl dust model is generally consistent with observations 
of Nova Cygni 1992, this conclusion would be overturned if the larger estimates of the reddening 
prove to be correct.^ Unfortunately, Nova Cygni 1992 is a less-than-ideal test of a dust model: as 
we have seen, there is uncertainty regarding the value of the reddening E{B — V) (i.e., the total 
amount of dust); in addition, there is uncertainty regarding the distribution of gas and dust along 
the sightline to Nova Cygni 1992 - it was plausible to consider that ~30% of the dust might be 
contributed by a cloud ~100pc from the nova. 

The ideal way to study scattering by interstellar dust is to employ an X-ray source well outside 
the galactic plane, so that most of the dust contributing to the scattering is at y < 0.1, in which case 
the halo angle 9 and scattering angle 65 are nearly the same for single-scattering. An extragalactic 
source would be ideal for this purpose. In this case, uncertainties regarding the precise location of 
the dust in the Galactic disk are unimportant (i.e., dust at y = and dust at y = 0.001 produce 
virtually identical halos) . We can hope that the great sensitivity of Chandra and XMM may make 
feasible observations of such out-of-plane sources. Alternatively, a source located in the galactic 
plane at 15° < |Z| < 75° would allow HI or CO observations, or optical absorption line spectroscopy, 
to locate the absorbing material using galactic rotation. 

The higher angular resolution of Chandra may allow observations at smaller halo angles, and 
the improved energy resolution will allow study of the energy spectrum of the scattered halo. 

8. Summary 

The principal results of this paper are as follows: 

1. We model the X-ray spectrum and light curve of Nova Cygni 1992 using the two-component 
model of Balman et al. (1998): optically-thin thermal plasma plus a 0-Ne enhanced white 



*X-rays scattered by dust very close to the nova would be lost in the point spread function. At 0.5 keV the median 
scattering angle Qm « lOOO", so most of the photons scattered by dust located at a; = 0.97 would be at halo angles 
6 < 30". Thus dust located at a; > 0.97 would not make an observable contribution to the X-ray halo. 
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dwarf atmosphere. 

2. We present the formalism for calculating X-ray scattering halos including multiple scattering. 
Time delay is treated exactly for single scattering and approximately for multiple scattering. 

3. We calculate the X-ray halo for Nova Cygni 1992, using the WDOl dust model, consisting 
of a mixture of amorphous silicate grains and carbonaceous/PAH grains. We compare the 
calculated X-ray halo with ROSAT PSPC imaging of the halo plus point source at 9 epochs. 

4. If the dust toward the nova is smoothly distributed (we consider an exponential density 
law as an example) we find that the WDOl dust model can quantitatively reproduce the 
observed halo intensity and angular profile provided the reddening E{B — V)'^ 0.20 mag, 
This is consistent with the value of E{B — V)k. 0.19 mag inferred from the observed Hq;/H/3 
emission line ratio at late times. 

5. If E{B — V)k, 0.20 mag, then the good agreement between the halo calculated for the WDOl 

grain model and observations of Nova Cygni 1992 contradicts previous claims that the dust 
toward Nova Cygni 1992 had to be either highly porous (Mathis et al. 1995) or include a 
substantial population of very large (o > 1/im) dust grains (Witt et al. 2001). 

6. Improved agreement between model and observation is obtained if ~ 30% of the dust is 
located in a cloud ~ 100 pc from the nova. 

7. The effects of time delays on the scattered halos depend on the distance to the source, which 
thus provides a method for distance determination to non-steady sources. The time delay 
of the halo relative to the nova is clearly visible at late times when the nova is in decline. 
The use of time-delay to determine the distance to Nova Cygni 1992 is discussed elsewhere 
(Draine & Tan 2003). 

8. It is hoped that future X-ray imaging by Chandra, XMM, or other telescopes will characterize 
the X-ray scattering halos around point sources where the dust column can be reliably esti- 
mated, and where there is information constraining the distribution of dust along the line of 
sight. Extragalactic X-ray point sources are ideal for this purpose. As we have shown here, 
such observations are capable of strongly testing models for interstellar dust. 
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